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1.0  FACTUAL  DATA 

1.1  Introduction 

The  study  supported  by  Contract  DA  36-039  SC  89099  is  intended  to  illuminate 
the  relationships  between  trc^ical  circulations  and  their  associated  distributions  of 
water  vapor  and  condensate.  A  reasonably  complete  description  of  these  relation¬ 
ships  must  consider  macro-  and  mesoscale  thermodynamic  influences  as  well  as 
cloud  physics  processes  and  the  kinematic  relationships  between  wind  and  water 
distributions. 

Although  there  has  been  substantial  progress  in  each  of  these  problem  areas, 
their  sjmthesis  has  been  considered  only  slightly.  The  present  effort  under  this 
contract  is  directed  toward  a  description  of  the  role  of  cloud  physics  processes, 
combined  with  modeled  tropical  circulations,  vapor  distributions,  and  condensation 
functions,  in  affecting  the  distributions  of  precipitation,  cloud,  and  vi^or. 

1.2  Collection  of  Unpublished  Results 
This  contract  has  enabled  us  to  finish  certain  investigations  that  have  been  in 

progress  for  the  past  two  years  and  to  prepare  the  findings  for  publication.  The 
typed  manuscript  was  submitted  to  the  Monthly  Weather  Review  on  July  31,  1962. 

Since  most  of  the  results  discussed  in  this  article  were  not  developed  under  this 
contract,  we  are  planning  not  to  print  it  as  a  scientific  report  under  this  contract. 

The  abstract  is  given  below. 

ELEMENTARY  THEORY  OF  ASSOCIATIONS 
BETWEEN  ATMOSPHERIC  MOTIONS 
AND  DISTRIBUTIONS  OF  WATER  CONTENT 

by 

Edwin  Kessler,  III 
ABSTRACT 

This  article  discusses  some  of  the  kinematic  relationships  between 
distributions  of  precipitation  and  air  motions,  based  on  analyses  of  conti¬ 
nuity  equations  embodying  the  principal  assumptions  that  (1)  water  v^or 
shares  the  motion  of  the  air  in  all  respects  and  (2)  condensate  shares 
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horizontal  air  motion,  but  falls  relative  to  air  at  a  speed  that  is  the 
same  for  all  the  particles  comprising  precipitation  at  a  particular 
time  and  height.  Among  topics  considered  are  the  steady-state  pre¬ 
cipitation  rate  in  a  saturated  horizontally  uniform  updraft  column;  the 
relationship  N  (V  w)  =  constant,  where  N  is  the  number  of  uniform 
particles  conq>rising  precipitation  at  a  particular  point  in  space  and 
time,  V  is  their  fall  speed,  and  w  is  the  updraft;  a  closed  formulation 
for  the  particle  concentration  and  size  in  precipitation  that  is  composed 
of  particles  of  the  same  size  at  any  point  in  space  and  time;  a  compar¬ 
ison  of  no  evaporation  and  instantaneous  ev^oration  of  condensate  in 
two  two-dimensional  model  wind  fields;  and  kinematic  descriptions  of 
the  radar -observed  generators  and  stalactites. 


1.3  Formal  Development  of  Solutions  of  the  One-dimensional  Continuity  Equation 
The  two  premises  noted  in  the  abstract  £d}ove  are  implicit  in  this  continuity 
equation  for  water  substance: 


8M  8M  8M  ,  8M  „  8V 

"ir  "  -  “  IT  -  '  TF  ■  "il  ■  “  ST  * 


(1) 


where  M  is  the  density  of  water  substance  in  all  its  phases  minus  the  saturation 

vapor  density;  G  =  -dQ  /dz,  where  Q  is  the  saturation  vapor  density;  V  is  the 

s  s 

fall  speed  relative  to  the  air  of  M  when  M  is  positive;  and  V  is  zero  when  M  is 
negative.  The  other  terms  have  their  usual  meanings. 

In  references  [3]  and  [4],  solutions  of  the  one-dimensional  form  of  Eq.  (1) 

for  an  inconqiressible  atmosphere  and  constant  fall  speed  of  condensate  are  given 

,  2  , 

for  the  limits  0  s  z  s  H  for  the  case  in  which  w  =  (4w  /H)  (z  -  z  /H).  The 

'  max  ' '  ' 

solutions  are  derived  from  two  total  differential  equations,  together  equivalent  to 
Eq.  (1),  which  describe  the  motion  and  development  of  individual  M-packets.  The 
solutions  presented  in  the  references  are  correct,  and  their  development  illuminates 
the  physical  processes  involved;  however,  it  seems  desirable  to  present  solutions 
developed  from  the  formal  logic  customarily  employed  by  mathematicians  for  such 
equations  and  within  the  historical  context  of  the  equations.  This  work  has  been 
done  by  Dr.  E.  A.  Newburg  and  will  be  presented  as  part  of  a  future  scientific  report. 
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1.4  A  Generalized  Model  that  IncludeB  Cloud  and  Rain 


In  the  earlier  studies  based  on  equations  like  (1).  cloud  is  omitted;  all  the 
condensed  water  is  assumed  to  exist  as  rain  or  snow  falling  at  velocity  V,  constant 
at  any  one  height. 

During  May,  June,  and  July,  1962,  we  have  generalized  the  earlier  theory  to 
include  an  elementary  modeling  of  the  bulk  effects  of  cloud  coalescence,  collection 
of  cloud  by  precipitation  and  evaporation  of  precipitation,  and  of  the  eStects  of  these 
processes  on  the  partitioning  of  water  substance  among  precipitation,  cloud,  and 
water  vapor.  The  following  are  the  principal  assumptions  inherent  in  our  preliminary 
mathematical  formulation. 

(1)  Cloud  is  condensed  water  that  has  no  ^preciable  fall  speed  relative  to 
the  air. 

(2)  Cloud  forms  in  rising  saturated  air,  and  ev£q>orates  in  saturated  descending 
air  at  the  rate  wG  =  -w  (dQ  /dz),  where  w  is  the  vertical  air  velocity  and  Q  is  the 

S  8 

saturation  v^or  density. 

(3)  Cloud  changes  to  raindrops  distributed  in  size  according  to  the  inverse 
exponential  distribution  of  Marshall  and  Palmer  [8]  at  the  rate  (m  -  a),  where 

and  a  are  constants  and  m  is  the  cloud's  water  content.  The  term  in  Eqs.  (2) 
and  (3)  that  embodies  this  assumption  is  called  the  autoconversion  term. 

(4)  Once  rain  is  formed,  it  combines  with  cloud  or  evaporates  in  subsaturated 
air,  in  accordance  with  approximate  forms  of  derived  physical  relationships.  The 
terms  in  Eqs.  (2)  and  (3)  that  accoimt  for  these  processes  are  the  accretion  of  cloud 
term  and  the  evaporation  term,  respectively. 

(5)  Although  an  inverse  exponential  distribution  of  precipitation  particles  is 
assumed  for  conq}uting  cloud  collection  and  evaporation,  the  vertical  mass  transport 
of  precipitation  is  based  on  the  terminal  fall  speed  of  the  drop  whose  diameter 
divides  the  distribution  into  parts  of  equal  water  content. 
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The  new  equations  for  vapor*  cloud,  and  precipitation  content  include  the 
macroscopic  effects  of  cloud  pl^sics  processes  and  incorporate  the  same  princ4>les 
of  continuity  discussed  in  the  references.  In  an  incompressible  atmosphere,  we 
have 


8m 

8t 


8M 

8t 


8m  8m  8m  , 

-  .  Kg„0.35^0.65^ 


8M  8M  ,  8M  „  8V 


-  a) 


+  Kj(m  -  a) 


(2) 


(3) 


where  m  may  be  positive,  negative,  or  zero  and  is  the  cloud  density  plus  the  vapor 
density  minus  the  saturation  vapor  density;  M  >  0  is  the  water  content  of  prec4>ita- 
tion;  E,  a  number  between  0  and  1,  assumed  independent  of  D,  is  the  efficiency  with 

which  precipitation  particles  collect  cloud;  K.  in  the  autoconversion  of  cloud  term 

-3  "1  ^  -4 

is  of  the  order  of  10  sec  i  Kg  ^  accretion  of  cloud  term  is  6.95  x  lO  when 

m  >  0  and  is  zero  otherwise;  K  in  the  evaporation  of  precipitation  term  is  1.35  x 
“6  ^ 

10  when  m  <  0  and  is  zero  otherwise;  and  V,  the  fall  speed  of  precipitation,  is 

,  T,  no  O  *“0.  125-  _0.125  y 

given  by  V  =  -OS.en^  M  m/sec. 

The  number  of  precipitation  particles  per  cubic  meter  of  air  per  unit  diameter 

interval  is  given  by 


-\D 

“  =  “o®  ' 

7  -4 

where  X  -  3.75/0^  (Atlas  (IJ).  The  order  of  magnitude  of  n^  is  10  m  . 


(4) 


1.5  Sample  Derivation 

As  an  exan^)le,  the  following  outlines  the  derivation  of  the  accretion  term  in 
Eqs.  (2)  and  (3). 

The  rate  at  which  volume  is  swept  out  by  one  drc^  of  diameter  D.  falling  at  a 

2  ^ 
speed  Vj  is  -irD^VjA  (where  V,  <  0),  and  the  rate  of  accumulation  of  cloud  water 
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by  a  single  precipitation  particle  is 


IT--—  '« 

where  is  the  efficiency  of  catch. 

The  rate  of  growth  of  the  li^d  water  content  in  the  entire  distribution  of 
precipitation  particles  is  given  by  integrating  Eq.  (5)  over  all  the  particles;  i.e., 

This  integration  is  acconq)li8hed  by  substituting  for  V  in  Eq.  (5)  according  to  the 
0.5  ^ 

relatibn  V  -  bO  ‘  .  due  to  Spilhaus  [9].  Then,  with  E  independent  of  D,  and  with  b 
sec  we  have 


4  TTEn^jin  /qD  e  dD, 


M  -  40  ^  r{3._5J  ... 

dt  4  j^3.5  • 

Equation  (9)  describes  the  collection  of  cloud  by  precipitation  and  is  obtained  from 
Eq.  (8)  by  substituting  for  \  its  equivalent  expression  in  terms  of  M: 

^  =  0.96  X  10  ^  EnQ’^^^mM®'®^®  gm  m  ®  sec  (9) 

The  derivations  of  the  fall-speed  equation  in  terms  of  n^  and  M  and  of  the 
evaporation  term  involve  integrals  similar  to  that  in  Eq.  (7).  The  ev£q)oration 
term  is  based  on  a  sinq)le  modeling  of  the  evaporation  functions  of  Kinzer  and 
Gunn  [6] . 

1.6  Method  of  Solution  of  Ecmatlons  of  the  Generalized  Model 

It  has  been  suggested  that  the  most  accurate  solution  of  Eq.  (1),  with  V  constant 
for  M  >  0,  is  obtained  by  considering  the  motion  and  development  of  M  along  its 
streamlines,  as  described  by  the  equivalent  total-differential  equations.  This  method 
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is  discussed  at  length  in  the  paper  whose  abstract  appears  in  Section  1.2.  Although 
such  Lagrangian  methods  of  solution  were  emphasized  in  the  original  proposal 
which  this  contract  now  supports,  these  methods  do  not  appear  feasible  in  connection 
with  the  solution  of  Eqs.  (2)  and  (3)  because  the  streanUines  are  not  fixed  as  in 
the  simpler  problem  and  application  of  the  method  would  involve  too  many  difficult 
and  lengthy  computations. 

We  have,  accordingly,  turned  again  to  finite-difference  methods  based  on  a 
rectangular  grid.  The  computer  program  prepared  by  Mr.  Gerald  Wickham  for  the 
one-dimensional  forms  of  Eqs.  (2)  and  (3)  follows  a  method  first  suggested  by 
P.  D.  Lax  [7] .  We  write  Eqs.  (2)  and  (3)  in  the  form 


8U  ^  §1  _  T, 

at  8z  ■ 

where  U,  F,  and  B  are  the  vectors 


LmJ  L  wm 


B(M)1 

LB(m)J’ 


(10) 


(11) 


and  the  B's  are  given  by 

dw 

B(M)  =  M  “  +  cloud  plqrsics  terms, 

oZ 

dw 

B(m)  =  m  —  +  wG  +  cloud  physics  terms. 

The  finite-difference  formulation  is 

U(z,  t+k)  =  U  (z,  t) 
av 


(12) 

(13) 


+  ^  {F(z-h,  t)  -  F(z  +  h,  t)  +  h{l-p)  [B(z  +  h)  +  B(z-h)]  +  2hpB(z)}, 

(14) 

where  k  is  the  time  increment,  h  is  the  height  increment,  p  is  an  e^qjerimental 
weighting  parameter  in  the  range  0  £  p  s  i,  and 

Uav(z.  t)  =  ^  lU(z+h,  t)  +  U(z-h,  t)J.  (15) 

At  a  lower  boundary,  we  use 

U(0,  t+k)  =  U(0,  t)  +  I  [F(0,  t)  -  F(h,  t)  +  hB(0,  t)],  (16) 


1 

1 

1 


6 


where  the  F'e  and  B's  do  not  include  terms  involving  w  or  8w/8z.  The  equations 
are  solved  under  the  mixed  conditions 

U  (z,  0)  =  g  (z)  (0  s  z  s  H). 

(17) 

M  (H.  t)  =  gg  (t)  (0  s  t  s  T). 

The  cloud  physics  terms  in  the  B's  have  been  modified  slightly  from  the  form 
given  in  Eqs.  (2)  and  (3),  so  that  their  first  derivatives  are  lipshitzian  [2].  The 
classical  Courant-Friedrichs-Lewy  stability  criterion  is  observed  in  choosing  the 
mesh;  the  domain  of  dependence  of  a  point  with  respect  to  the  differential  equation 
always  lies  within  the  domain  of  dependence  of  this  point  with  respect  to  the  differ¬ 
ence  equation. 

1.7  Sample  Solution 

First  problems  solved  by  the  IBM  7090  conq>uter  included  those  simple  one- 
dimensional  cases  discussed  in  the  references,  whose  solutions  are  accurately  known. 
In  this  progress  report,  we  note  merely  that  the  computer  answers  are  highly  gratify¬ 
ing;  the  finite-difference  scheme  is  distinctly  superior  to  that  used  in  the  earlier 
work. 

Figure  1  shows  some  of  the  results  obtained  from  the  first  run  of  the  computer 
on  a  problem  that  included  conversion  of  cloud  to  precipitation  and  collection  of 
cloud  by  precipitation.  The  solution  demonstrates  many  features  of  a  shower  which 
appear  analogous  to  events  in  the  real  atmosphere,  although  the  input  parameters  for 
this  first  trial  constituted  a  relatively  incon4)lete  simulation  of  the  real  tropical 
atmosphere.  Further  discussion  of  this  and  other  solutions  will  be  included  in  the 
next  progress  report. 
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2.0  CONCLUSIONS 


The  presently  available  one-dimensional  computer  program  should  contribute 
greatly  to  evaluation  of  the  relationships  between  the  developments  of  precipitation 
and  the  strengths  of  condensation,  coalescence,  accretion,  and  evaporation  processes. 
The  analysis  of  one-dimensional  solutions  should  be  vigorously  pursued. 

The  study  of  budget  parameters  will  require  computer  programs  to  solve  the 
continuity  equations  in  two  and  three  dimensions  [3,  4]. 
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3.0  PROGRAM  FOR  NEXT  INTERVAL 

Theoretical  and  numerical  approaches  will  be  used  to  study  solutions  of  the 
one-dimensional  equations.  Among  the  mai^  questions  to  be  addressed  are: 

(1)  Under  what  conditions,  If  az^,  are  steady  updrafts  and  cloud  pl^slcs 
parameters  associated  with  oscillatory  behavior  of  the  cloud  and  precipitation 
distributions? 

(2)  How  is  the  time  of  onset  of  precipitation  related  to  the  strengths  of  the 
coalescence  and  accretion  processes?  How  do  variations  of  these  processes  affect 
the  peak  rate  of  rainfall  and  its  rise  time? 

(3)  How  do  the  rate  and  duration  of  surface  rainfall  vary  with  the  strength  and 
duration  of  the  iq)drafts  ? 

(4)  How  do  the  several  variables  relate  to  the  persistence  of  cloud  aloft? 

(5)  Can  the  shape  and  intensify  of  the  wind  field  be  deduced  from  observations 
of  the  cloud  and  precipitation  distributions? 
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